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This paper deals with adaptive transmitting arrays in which the excita- 
tions of the elements are varied in response to a pilot field incident on the 
array from a distant source. General theorems, some quite simple, are ob- 
tained relating to optimal power transfer from an adaptive array in an arbi- 
trary reciprocal medium to either a single receiver or a receiving array. We 
assume first that (he amplitudes and phases of the transmitting elements are 
separately adjustable, and afterward that only the phases are adjustable. 
The results involve in particular the matrix which represents the pilot fields 
produced at the elements of the transmitting array by currents at the locations 
of the receiving elements. In some important special cases, optimal power 
transfer results from making the phase of each transmitting element equal to 
the negative of the phase of the pilot field at that element. 

We also consider the dynamic behavior of two adaptive arrays which simul- 
taneously transmit and receive, the phases on transmission being made equal 
to the negatives of the received phases. Analysis of an idealized model indi- 
cates that the arrays will reach a unique steady state which is in practical 
cases identical with or very close to the condition for optimal power transfer. 
Some numerical simulations of 2- and 8-element interacting arrays have 
been made to show how such arrays approach an essentially steady state 
under moderately realistic assumptions. 

I. INTRODUCTION AND SUMMARY 

A number of recent papers 1 ■"- have dealt with adaptive antenna arrays, 
also called self-steering or retrodirective arrays. In an adaptive trans- 
mitting array, the excitations of the individual elements are electroni- 
cally varied in response to a pilot field incident on the array from a 
distant terminal, in order to steer the beam to the terminal which is origin- 
ating the pilot signal. It is easy to see that in free space the required 
steering can be accomplished by making the phase of each element equal 
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to the negative of the phase of the pilot beam at the given element. 
Cutler, Kompfner, and Tillotson 1 and others 2 have shown how phase 
reversal can be obtained using frequency conversion techniques. 

This paper deals with general adaptive transmitting arrays in an 
arbitrary environment. The transmission medium need not be homogene- 
ous or isotropic, but it is assumed to be linear and symmetric and to be 
time-invariant, at least over intervals comparable to the propagation 
time between transmitter and receiver. We are concerned in particular 
with the conditions for optimal power transfer from an adaptive trans- 
mitting array to either a single receiver or a receiving array. We shall also 
investigate the transient and steady-state behavior of two interacting 
adaptive arrays, each of which simultaneously transmits and adjusts 
the excitations of its elements in response to the field received from the 
other array. 

In Section II we consider an adaptive array in which the amplitudes 
and phases of the excitations of the individual elements are separately 
adjustable, but the total radiated power is fixed. Such an array is easier 
to treat mathematically than one in which the excitation amplitudes are 
all fixed and only the phases are variable, even though the latter array 
might be easier to build. In the most general case, the power radiated by 
the transmitting array is a positive definite Hermitian form in the ele- 
ment excitations, and the received power is a positive definite Hermitian 
form in the electric fields at the elements of the receiving array. The dis- 
tribution of excitations which maximizes the ratio of received to radiated 
power is the eigenvector corresponding to the largest eigenvalue of a 
certain pencil of Hermitian matrices. The matrices in question are con- 
structed from the impedance matrix of the transmitting array, the ad- 
mittance matrix of the receiving array, and a Green's function matrix 
of pilot fields produced at the transmitting elements by currents at the 
receiver locations. The results simplify considerably if the elements of 
each array are uncoupled and are identical among themselves. In par- 
ticular, if the receiver consists of but a single element, and if the trans- 
mitter elements are identical and uncoupled, then the optimal excitation 
of each element is merely proportional to the complex conjugate of the 
pilot field at that element. 

Section III contains a brief discussion of the problem of maximizing 
the power transferred from an arbitrary transmitting array to an arbi- 
trary receiving array when the excitation amplitudes are fixed and only 
the phases are adjustable. Maximum power is always conveyed to a 
single receiver by reversing the phase of the pilot field at each element of 
the transmitting array. This phase reversal principle, first recognized for 
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free-space transmission, is thus shown to be valid for an arbitrary trans- 
mission medium. For multielement receivers an explicit solution is not 
given ; but an example shows that even when each array consists of iden- 
tical, uncoupled elements, maximum power transfer does not generally 
correspond to reversing the phase of the total pilot field at the trans- 
mitter elements. 

Section IV deals with the interaction of two adaptive arrays or, in 
principle, the interaction of an adaptive radar array with itself. A mathe- 
matical model is set up, in which each array transmits constant power and 
continuously adjusts the excitations of its own elements to be propor- 
tional to the complex conjugate of the incident field. A single delay time 
is taken to represent the transmission delay between the two arrays. The 
transient behavior of this model turns out to be quite simple, and it is 
shown that in general, excluding mathematically pathological cases, the 
two arrays reach an equilibrium configuration which depends only on the 
Green's function (pilot field) matrix corresponding to the given geome- 
try and transmission medium. In the most general case the equilibrium 
configuration is not the same as the condition for optimal power transfer 
derived in Section II ; but it is the condition for optimal power transfer 
in the important special case when the elements of each array are identi- 
cal among themselves and the interelement coupling is zero. If the ele- 
ments are nearly identical and the mutual impedances are small com- 
pared to the self-impedances, then the equilibrium configuration should 
be nearly the same as the configuration for optimal power transfer. 

Numerical simulations of the transient behavior of 2- and 3 -element 
interacting adaptive arrays are described in Section V, both for the case 
of simultaneous phase and amplitude variations, and for the case of phase 
variations only. The simulations also include the effects of small dif- 
ferences in the interelement delay times compared to the average delay 
between the arrays. Random choices are made for the elements of the 
Green's function matrix and for all pairs of interelement delays. Simula- 
tions of 50 pairs of 2 -element arrays and 25 pairs of 3-element arrays 
indicate that arrays with only phase adjustment approach a steady 
state about as quickly as arrays with both phase and amplitude adjust- 
ment (of course, the two steady states are not the same). Interelement 
delay differences which are small compared to the average interelement 
delay produce small fluctuations about the steady state which would be 
achieved for equal delays. 

The results obtained in this paper depend only on the linearity, sym- 
metry, and time-invariance of the transmission medium; in particular, 
they do not involve calculating any antenna patterns. Pattern calcula- 
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tions would be necessary if one wished to get numerical values for maxi- 
mal power transfer, or to estimate the radiated fields in unwanted direc- 
tions. Furthermore, the analysis is essentially for a single frequency; 
variations in phase and amplitude are assumed to be very slow compared 
to the transmission times involved. It would be worthwhile to study the 
behavior of adaptive arrays over a finite frequency band, but such a 
study is outside the scope of the present paper. 

II. OPTIMAL POWER TRANSFER BETWEEN ARBITRARY ANTENNA ARRAYS 

Consider a transmitting array and a receiving array embedded in an 
arbitrary linear, time-invariant medium, as in Fig. 1. The medium maj r 
be inhomogeneous and anisotropic, but the permeability, permittivity, 
and conductivity tensors at any point are assumed to be symmetric. 
(This rules out ferrites and plasmas in the presence of a magnetic field.) 
All fields are assumed to be time-harmonic with angular frequency co, 
the time dependence exp iwt being suppressed. For simplicity the in- 
dividual radiators and receivers are taken to be elemental electric dipoles, 
although they could equally well be elemental current loops. The assump- 
tion of dipole sources is not a major restriction, since the dipoles could 
be used, for example, together with microwave circuitry to feed aperture- 
type radiators such as elemental horns. 

Let the transmitting array have M elements and let the complex 
excitation of the ith element be h,i . Physically h,i may be regarded 
as the electric moment of an elemental current, having the dimensions 
of ampere-meters. The M-component vector 

II = (/l.l , /l.2 , ••• ,/..«), (1) 

whose components are complex scalars, will be called the excitation of 
Array 1. Similarly let the receiving array have JV elements, and let the 
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Fig. 1 — Schematic representation of arbitrary transmitting and receiving 
arrays of electric dipoles. 
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complex electric field component at the location and in the direction of 
the jth element be E 2 j . Then the electric field at Array 2 is the ^-com- 
ponent vector 

E 2 = (E 2A ,E 2 , 2 , ••• ,E 2 , N ). (2) 

If Array 2 is transmitting and Array 1 is receiving, we define the vectors 
I 2 and E] in an analogous way. 

The total power P T radiated by the transmitter over all space is given 
by the Hermitian form 

F T = HZiIi,Ii), (3) 

where Z x is an M X M positive definite Hermitian matrix and (x,y) 
represents the scalar product of two vectors x and y. A fuller discussion 
of notation and of the properties of Hermitian forms is given in Appendix 
A. 

In principle the radiation impedance matrix Zi may be determined 
from experimental measurements, or it may be calculated from the 
fields of the radiating elements. For example, if the field due to unit 
excitation of the zth element at a great distance R from all currents and 
material media is &i ti /R, then by integrating the Poynting vector due 
to the whole array over a large sphere we find that the total radiated 
power is given by an expression of the form (3), with 

Zx.it = - f &i,i-&i,j*dtt, i = 1,2, • • -, M; j = 1,2, • • •, M, (4) 

where r\ is the characteristic impedance of free space and dtt is an element 
of solid angle. 

We now assume that the power P R received by Array 2 is given by a 
Hermitian form in Eg , the electric field which would exist at Array 2 if 
its elements were open-circuited. Thus we write 

P R = i(Y 2 E,,E 2 ), (5) 

where Y 2 is an N X iV positive definite Hermitian matrix. Equation (5) 
is equivalent to the assumption that the transmitter field is independent 
of whether or not currents are flowing in the elements of Array 2, i.e., 
that the back reaction of Array 2 on Array 1 is negligible. This will be a 
very good approximation in the practical case where the arrays are far 
apart, so that P R is a very small fraction of P T ■ 

The field at Array 2 is related to the excitation of Array 1 by the Green's 
function matrix T; thus 

E 2 = rL , (6) 
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where r is an N X M matrix and r,-> represents the field at the ith ele- 
ment of Array 2 due to unit excitation of the jth. element of Array 1. A 
basic reciprocity theorem for linear, time-invariant, symmetric media, 
proved in Appendix B, guarantees that r u also represents the field at 
the jth element of Array 1 due to unit excitation of the ith element of 
Array 2. Thus if Array 2 has the excitation I 2 , the field at Array 1 is 
given by 

E, = r'l 2 , (7) 

where r' is the transpose of r. 

From (5) and (6), the power received by Array 2 is 

P R - |(Y 2 rli , no = |(i*T,l% , 10, (8) 

where r + is the adjoint ( = conjugate transpose) of r. We wish to maxi- 
mize the ratio of received power to transmitted power, which is 

i\ = (rnwyy (9) 

Pt ' {Zxlik) 

But the right side is the quotient of two Hermitian forms in which the 
denominator is positive definite; and it is well known (see Appendix A) 
that the maximum value of the quotient is the largest eigenvalue X M of 
the pencil of matrices r f Y 2 r — XZi . The desired eigenvalue is the largest 
root of the equation 

det (r + Y 2 r - XZi) = 0; (10) 

and the corresponding eigenvector, which maximizes the right side of 
(9), is any nonzero solution of the system of equations 

r^rli - XjiZiIi = o. (li) 

The foregoing equations simplify in a special case which will be im- 
portant in what follows, namely when all the self-impedances and self- 
admittances are equal and all the mutual impedances and admittances 
are zero. In this case we may write 

Zi^fijljf, Y 2 = C 2 1„, (12) 

where Ri and (7 2 are real scalars and 1m and 1a- are unit matrices of orders 
M and N respectively. Then (9) becomes 

P R _ g 2 (rWjk) (13) 

Pt Ri (Ii,Ii) ' 
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and the maximum value of the ratio is proportional to the largest eigen- 
value of the matrix r + r, that is, the largest root Aj, of 

det (r + r - XI m) = 0. (14) 

The excitation corresponding to maximum power transfer is any non- 
zero solution of 

r f rii-AjJi = o. (is) 

The optimal transmitter excitation given by (14) and (15) is one 
which can exist when both arrays are transmitting and the excitation of 
each element of each array is proportional to the complex conjugate of 
the field incident on the element from the other array. Suppose, for ex- 
ample, that 

I, = J/jE,*, I- = M&\ (16) 

where Mi and M 2 are complex scalars and Ei and E 2 are related to I 2 
and Ii by the Green's function matrix, as in (6) and (7). Then it is easy 
to show that Ii must satisfy 

I t = M 1 M 2 *r t rl 1 . (17) 

A nonvanishing solution of (17) exists if and only if 

MiMt* = 1/X, (18) 

where X is an eigenvalue of r f r, that is, a root of the determinantal equa- 
tion (14). Although steady-state excitations satisfying (16) are mathe- 
matically possible when X is any eigenvalue of r + r, it is shown in Sec- 
tion IV that the system is unstable unless X is the largest eigenvalue, and 
that the excitations corresponding to the largest eigenvalue are in fact 
the excitations toward which two interacting adaptive arrays tend. 

We shall now consider optimal power transfer in the special case where 
Array 2 consists of but a single receiver. Then Z x is an M X M matrix, 
Y 2 is a 1 X 1 matrix, i.e., a scalar quantity G% , and r is a 1 X M matrix. 
Equation (10) is therefore equivalent to 

det (xrAf'z, - r f r) = det (\G.r%) det (1,, - [xGfXr 1 r*r) 

= X A/ &T , 'detZ 1 det (1, - r{\G 2 ~%r 1 t [ ) 

(19) 
= (X/G' 2 )- U_1 (X/G 2 - rzr'r*) det Zx 

= 0, 

where in the second step we have made use of a lemma due to Sandberg,' 
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which is stated and proved in Appendix C. It follows that the only non- 
zero eigenvalue of (19) is 

\ M = G 2 rzr l r\ (20) 

The corresponding eigenvector is 

i! (J/) = zrV (2i) 

up to a constant factor, since it is easy to see that 

i*ftHk <J0 - \*Zili (J0 = r^rzrV - x^zr'r* 

. + (22) 

= r*Xj, - \ M r f = o. 

If the elements of the transmitting array are uncoupled, then the 
mutual radiation impedances vanish, and Zi and Z{~ 1 are diagonal mat- 
rices. The optimal excitation of the jth element is then 

I 1J = Tif/Bu, 3 = 1,2,--- ,M, (23) 

up to a constant factor, where Ti* is the complex conjugate of the pilot 
field produced at the jth element by a dipole at the location of the re- 
ceiver, and Rjj is the radiation resistance of the jth element. If all the 
radiation resistances are equal, then since the eigenvector is determined 
only up to a multiplicative constant, we may take 

/i.i-r w * f j= 1,2, ••• ,M; (24) 

in other words, the excitation is merely proportional to the complex 
conjugate of the pilot field. 

We have shown that if the mutual radiation impedances of the trans- 
mitter elements are zero, then the field at the receiver is maximized, for 
constant radiated power, when the phase of the excitation of each trans- 
mitter element is the negative of the phase of the pilot field. If, however, 
the transmitter elements are coupled by their radiation fields, so that the 
impedance matrix Zi is not diagonal, the optimal excitations are given 
by (21) and do not generally satisfy the phase reversal condition. 

III. OPTIMAL POWER TRANSFER WITH PHASE ADJUSTMENTS ONLY 

If the amplitudes of the transmitter excitations are fixed but the phases 
are adjustable, we wish to maximize the received power, 

p R = Kr^ri, ,i,) = j(T,Hk , no, (25) 

when 
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h.i-r/"*, j=l,2,--,M, (26) 

and the rj arc fixed but the 0> are at our disposal. 

If there is only one receiver (N = 1), the solution is immediate. We 
have to maximize 



Pr = ¥h 



E iVw 



j=i 



(27) 



by adjusting the phases of the l\,j ; and it is clear that the modulus of the 
sum will be greatest when the phases of all the summands are equal, that 
is, when 

arg h,i = Bj = -arg Tu + constant, ; = 1, 2, • • • , M. (28) 

In other words, the phase of the jth transmitting element should be the 
negative of the phase of the pilot field produced at that element by a 
radiating clement at the position of the receiver. This result is independ- 
ent of the nature of the transmission medium, subject only to the re- 
quirements of linearity, time-invariance, and symmetry, and it is in- 
dependent of the position of the receiver relative to the transmitting 
array. 

For a two-element transmitter (M = 2) and an arbitrary receiving 
array, we have 

2P I{ = (I*Y,rW + (r t Y 2 r) 12 r 1 r 2 e" 92 - 9j5 

+ ( r t Y 2 r) 2l r 1 r 2 e ,(fll ^ 2) + ( r t Y 2 r) 22 r 2 2 . 

Since r + Y 2 r is Hermitian, the right side of (29) is maximized by taking 

arg /,., - arg /•>,, = 0, - e, = arg (r + Y 2 r) 12 , (30) 

and this is the condition for optimal power transfer if the transmitting 
array has only two elements. 

A complete analytic solution of the problem of maximizing P R for an 
arbitrary transmitter with il/ ^ '.i and an arbitrary receiver with N ^ 2 
has not been found, although since P lt is a continuous, periodic function 
of each of the Oj , it is obvious that a maximum exists and could be lo- 
cated as accurately as desired by an iterative numerical procedure. 

In contrast to the situation for arrays with both amplitudes and phases 
adjustable, the condition for optimal power transfer between multiele- 
ment arrays with fixed excitation amplitudes is not generally satisfied 
by making the phase of each element equal to the negative of the phase 
of the field incident from the other array, even if it be assumed that the 
array elements are uncoupled and are identical among themselves. As a 
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counterexample, consider the case in which each array has two elements, 
and Y 2 is a multiple of the unit matrix. Let the element currents be 

h,\ = re \ i 2 ,i = pe r, 

(31) 
/i, 2 = re te \ h.2 = pe 9 \ 

where r and p are real and positive, and the phases are at our disposal. 
Equations (6) and (7) give 

# 2 ,i = r(r 11 e"' 1 + iV 2 ), 

En = r(r 21 e' 91 + rue**), 

(32) 

E ltl = p(r u e"" + r,ic ,M ), 

E lti = p(r 12 e ,Vl + r 22 e ,V2 ). 

The phase reversal condition leads to the following pair of simultaneous 
equations : 



d, - 6-2 - arg «i.i - arg E UI = arg p <(yi - w) + ^ 

„ r n e < (,1 "* ,) + r 22 

Vi - «0a = arg #2,2 - arg A 2 ,i = arg i(g e) 

I ne ■ -j- 1 12 



(33; 



On the other hand, the condition (30) for maximum power transfer re- 
duces to 

e t - 02 = arg (r^r)u = arg (ru*r M + r 2 i*r 22 ). (34) 

Since r is an essentially arbitrary complex matrix, equations (33) are 
not equivalent to (34), although it is possible that in practical cases the 
two conditions will yield values of 0i - d-> which do not differ by very 
much. 

IV. DYNAMIC BEHAVIOR OF INTERACTING ADAPTIVE ARRAYS 

In this section we set up a simple model of the dynamic behavior of 
two adaptive arrays, each of which continuously adjusts the excitations 
of its own elements in response to the fields from the other array. In 
principle the same equations would apply to a single array interacting 
with itself, as a combined radar transmitter and receiver. The funda- 
mental assumption is that the amplitudes and phases of the element cur- 
rents vary so slowly, compared with the transmission time between the 
arrays, that a single-frequency analysis is valid. 

Since in this model the excitation of an adaptive array is indeterminate 
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in the absence of an external field, we have to use an auxiliary antenna 
or beacon to turn the system on. The steps are as follows: First the bea- 
con is turned on, illuminating at least Array 1. Then Array 1 is turned on, 
Array 2 is turned on, and the beacon is turned off, leaving Arrays 1 and 2 to 
interact only with each other. It is convenient to assume that when a 
transmitter is switched on or off, its radiated power changes continu- 
ously, during a finite time interval, from one steady-state value to an- 
other. 

First we consider arrays in which the excitation of each element is pro- 
portional to the complex conjugate of the field incident on that element, 
and the total radiated power is a prescribed function of time. Thus let 
Ii(/) and I 2 (0 be the (slowly varying) complex excitations of the two 
arrays. We assume that the dynamic behavior of the arrays is described 
by the following equations : 

h.j(t) = Mi(0e Wl [&,,*/..**(/ - r kj - n) + &./■« - n)J, 

j = 1,2, ••• ,M; 

hAO = mi*)*" 9 [g T>*li/(t ~ *« - «) + **.**« " *) ]« (36) 

k = 1, 2, • • • , N. 

In these equations, Bi(0 and B 2 (/) are the beacon fields, if any, at Arrays 
1 and 2, r*,- is the transmission delay between thejth element of Array 1 
and the fcth element of Array 2, u and t 2 are constant time delays in the 
amplifiers of Arrays 1 and 2, t?i and d 2 are constant phase shifts, and 
Hi(t) and no(t) are real normalization factors determined by 

h(ZMt)Mt)) = Pi(t), (37) 

h(ZMt)MD) = Pi(0, (38) 

where the radiated powers /MO and P 2 (t) are given functions of time. 
Similarly, the equations describing two arrays in which the excitation 
amplitudes | Iij(t) \ and | / Sl *(0 | are prescribed functions of time, while 
the phases are continuously adjusted to satisfy the phase reversal condi- 
tion, are as follows: 

arg hj(t) = 0i - arg I fj V kj hAt - r ki - n) + B llS (t - n)J, 

j = 1, 2, • • • , M; 
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arg h,k(t) = #2 - arg X V k jhM ~ thj - r%) + B 2lk (t - t 2 ) \, 

Li=i J (40) 

fc - 1, 2, ■ • • , N. 

Before undertaking numerical simulations of the dynamic behavior 
of adaptive arrays, we consider an example which can be handled ana- 
lytically, namely the special case of two power-limited arrays in which 
all the interelement delay times are equal. We obtain this case from (35) 
and (36) by setting 

Tkj = t 3 = constant. (41) 

If we assume for simplicity that the beacon has been turned off, (35) 
and (36) take the form 

h(t) = m(0e* l rV(< - n - n), (42) 

Ut) = ^(t)e^-r%*(t - r, - r 2 ). (43) 

Eliminating I 2 yields 

1,(0 = »(0«*I*Iti(« - r), (44) 

where 

= 0! - 2 , T = Tl + T , + 2r 3 , (45) 

and the normalizing factor ni(t) may be expressed, if needed, in terms 
of the radiated power Pi(t) by (37). 

The Hermitian matrix r f r is at least positive semidefinite and will 
have M real eigenvalues. We suppose that the eigenvalues are numbered 
in order of increasing size and that the largest eigenvalue is unique; 
that is, 

^ Xi ^ X 2 ;g ••• ^ Xa,_, < \ M . (46) 

The corresponding eigenvectors z (1) , z (2) , ■ • ■ , z {M) satisfy 

r f rz (0 = X,-z (0 , i - 1, 2, ... ,M, (47) 

and may be taken as ortho normal, i.e., 

(z (i) ,z U) ) =««,. (48) 

Let Ii(t) be expanded in terms of the z (,) , with coefficients depending, 
of course, on time; thus 

li(0 = T,Ci(t)z w . (49) 
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Repeated application of (44) and (47) gives 

h(t + nr) = N n (t)e in * EcXOA/'z 01 , (50) 

3=1 

where N„(t) is again a normalization factor chosen to satisfy (37). 

Now suppose that there is an interval of length t, say to ^ t < t + r, 
in which all the c,(/) are hounded and c.u(t) is bounded away from zero. 
It follows from (46) that the term in X. w " in (50) will eventually dominate 
all the others, and we shall have 

_ . N [2Pi(0] J exp [in& + arg c M (t)] {M) 

Ut + n7) B ?. ~ (Z^/"')' ' (51) 

for t ^ t < t + r. 

It is easy to verify that the phase of li(0, as given by (51), is continuous 
at t = to + (n + 1)t if the phase of Cat(0 is continuous at t = t + t. 
We have just proved that two power-limited adaptive arrays with 
equal interelement delays will reach an equilibrium state in which the 
excitation of Array 1 is proportional to the eigenvector belonging to the 
largest eigenvalue of r + r. Similarly, the equilibrium excitation of Array 
2 is proportional to the eigenvector belonging to the largest eigenvalue of 
rr + (again X M )- But it was shown in Section II that the eigenvector be- 
longing to the largest eigenvalue of r f r corresponds to maximum power 
transfer whenZi and Y 2 are multiples of the unit matrix; that is, when 
all the mutual impedances and admittances are zero and all the self- 
impedances and self -admittances are equal. If this condition is approxi- 
mately satisfied, as will often be the case in practice, then the equilibrium 
excitation should be nearly the same as the excitation for optimal power 

transfer. 

We observe that the equilibrium excitation is unique except in the 
pathological case where r + r has two or more equal eigenvalues which are 
larger than all the rest, Steady states in which the current distribution 
corresponds to one of the smaller eigenvalues of r f r are mathematically 
possible, but are unstable. Also, if the arrays are moved with respect to 
each other or if the transmission medium changes (either case would 
correspond to changing the Green's function matrix), the final equilib- 
rium state depends only on the final positions of the arrays and the prop- 
erties of the transmission medium, and not at all on how the situation 

was reached. 

It should be pointed out that the foregoing argument does not apply, 
at least in its present form, to fixed-amplitude arrays with only phases 
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adjustable. Clearly there are extreme cases in which the behavior of 
fixed-amplitude arrays will be qualitatively different from that of power- 
limited arrays. For example, if r is diagonal, so that each element of Ar- 
ray 1 is coupled to only one element of Array 2, then power-limited arrays 
will ultimately cut out all of the elements except for the pair which is 
most closely coupled ; but the elements of fixed-amplitude arrays will go 
on indefinitely talking to each other in pairs, with no particular phase 
relationship between the elements of different pairs. Nevertheless, the 
numerical simulations of the next section indicate that in typical cases 
fixed-amplitude arrays do settle down to a steady state about as quickly 
as power-limited ones. As yet, however, no mathematical theorem has 
been proved about the steady-state behavior of fixed-amplitude arrays. 

V. NUMERICAL SIMULATIONS 

Because the equations describing the dynamic behavior of interacting 
adaptive arrays generally do not lend themselves to analytic treatment, 
we have made a few numerical simulations of 2- and 3-element arrays on 
an IBM 7094, in order to get some feeling for the possible behavior of 
interacting adaptive arrays in practice. Since these simulations were 
only computational experiments, no physical significance is to be at- 
tached to the specific numerical results. 

We shall first describe the method of simulation, then show the out- 
come of a typical calculation, and finally summarize the results of the 
whole study. 

For each simulation we selected the elements of the 2 X 2 or 3 X 3 
matrix r according to the following scheme : We set 

IV = Gjtexp (-iirrijk/h), (52) 

where Gry* was a random number selected with equal probability from 
the set { 1,2,4,8} , and n# was selected with equal probability from the 
set { — 5, —4, • • • , 5} . For the interelement delay times we took 

Tj , = 20 4- n ik . (53) 

As a matter of interest, we also computed the Hermitian matrix r + r 
and its eigenvalues and eigenvectors. 

If the time delays are all commensurable, (35) and (36) or (39) and 
(40) can easily be solved recursively on a digital computer. To start the 
system off, Array 1 was supposed to be illuminated initially by a con- 
stant beacon field. Array 2 was turned on linearly during a period of 20 
time units, and the beacon was turned off linearly during a similar period. 
Four cases were run with each choice of T. 



ADAPTIVE ANTENNA ARRAYS 37 

Case I. Power limited, equal delays. The condition 

(I,,Ii) = (I 2 ,I 2 ) = 1 (54) 

was imposed, and all interelement delays were set equal to 20 units. This 
case must approach a steady state, according to Section IV, provided 
only that the largest eigenvalue of r f r is unique. 

Case II. Power limited, unequal delays. Same as Case I, except that 
the time delays given by (53) were used. 

Case III. Fixed amplitudes, equal delays. The condition 

|j 1§y | = |/ 2 ,, | = 1 (55) 

was imposed, and all interelement delays were set equal to 20 units. 

Case IV. Fixed amplitudes, unequal delays. Same as Case III, except 
that the time delays given by (53) were used. 

In a typical run, the random number generator produced : 

/-I -1 -2\ 
(»») = -5 -5 , (56) 

\-l 1 3/ 

/8/36 2/36° 8 /72° \ 
0»=( 2 ZQ! 8 /180° 1 /180° . (57) 

\8/36° 8 / -36° 2/ -108 °/ 

It follows that 

/132.0/0 64.Q /-72.Q 46.4 /37.5 

r t r = ( 64.0 772.0° 132.0/0° 26.5 /-12.7 

\ 46.4 /-37.5° 26.5 /12.7° 69.0/0 

X, = 20.6, X 2 = 109.8, Xa = 202.6. 
The eigenvector corresponding to X 3 is 

z (3) = (0.719/0°, 0.656 /64.5° , 0.22 9/ -6. 2° ). (60) 

Figs. 2 through 5 show the results of running Cases I through IV over 
the time interval ^ / ^ 400. In the figures the phases are referred to 
the phase of the first element of each array, and the following notation 
is vised: 

i! = (A 1 ,A 2 /^2, A 3 M), (61) 

I 2 = (B 2 ,B,/l2,BM. (62) 
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The initial behavior of the two arrays depends on the particular way 
in which they were turned on and is of no great importance; what we are 
really interested in is the behavior at large times. In Cases I and III 
(Figs. 2 and 4), under the assumption of equal interelement delay times, 
the system appears to settle down to a perfectly steady state. It is easy 
to verify that in Case I the steady-state excitation of Array 1 corresponds 
to the eigenvector z <3) given by (60). On the other hand, in Cases II and 
IV (Figs. 3 and 5), where the delays are not all equal, the array excita- 
tions continue to show small residual fluctuations about the steady-state 
solutions of Cases I and III. These fluctuations are quite apparent in the 
original plots from which the present figures were redrawn. 

In the numerical study, 50 pairs of 2-element arrays were simulated 
and four cases run for each pair. The ratio of eigenvalues X2/X1 of r + r 
ranged from 293.5 to 1.385. As expected, the larger values of X 2 /Xi gener- 
ally produced quicker convergence; but only one case, out of all those 
tried, failed to reach essentially steady values by t = 400. In this par- 
ticular example X2/X1 was 7.37, and the interelement delays happened to 
range all the way from 15 to 25. Cases I, II, and III settled down rela- 
tively quickly, but Case IV (fixed amplitudes, unequal delays) went into 
a large-amplitude oscillation which was obviously not dying out at 
t = 1000. A similar, subsequent run in which the extreme interelement 
delay times were changed to 16 and 24 settled down normally. 

Twenty-five pairs of 3-element arrays were simulated, with eigenvalue 
ratios X 3 /X 2 ranging from 26.02 to 1.458. Every one of these cases ap- 
peared to have reached an essentially steady state at t = 400. The ex- 
ample shown in Figs. 2 through 5 is entirely typical. 

From the numerical simulations it is clear that sufficiently large delay 
differences (perhaps ±25 per cent of the average delay time) can make 
a pair of interacting adaptive arrays fail to settle down. We conjecture, 
however, that the arrays will always reach an essentially steady state if 
the delay differences are a sufficiently small fraction of the average delay. 
Conceivably one could put bounds on the fluctuations as a function of 
the deviations of the delays from the mean delay, but a more practical 
approach might be to do some experiments with real adaptive arrays. 
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Fig. 2 — Transient behavior of amplitudes and phases in power-limited adap- 
tive arrays with equal interelement delays: (a) Array 1, (b) Array 2. 
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Fig. 3 — Transient behavior of amplitudes and phases in power-limited adap- 
tive arrays with unequal interelement delays: (a) Array 1, (b) Array 2. 
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Fig. 4 — Transient behavior of phases in fixed-amplitude adaptive arrays 
with equal interelement delays. 



APPENDIX A 



Vectors, Matrices, and Hermitian Forms 

We summarize here the notation used in this paper, as well as some 
properties of Hermitian forms which are proved in textbooks like that 
of Gantmacher. 4 

A vector in n-diniensional complex space is an ordered array of n 
complex numbers: 

X = (X\ ,X 2 , • ■■ ,Xn). (63) 

The scalar product of the vectors x and y is written (x,y) and is defined 

by 



(x,y) = £.Ct2/.*, 
where the asterisk denotes complex conjugate. 



(64) 
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Fig. 5 — Transient behavior of phases in fixed-amplitude adaptive arrays 
with unequal interelement delays. 

A matrix is an m X n array of complex numbers : 

A = (A u ), i=l, 2, ...,m; j - 1, 2, • • • , n. (65) 

Associated with a given matrix are the following matrices: 

Conjugate (A ),-; = A i} *, 

Transpose (A')»j = A >,- , (66) 

Adjoint (A f )ij = Aj*. 

Note that this definition of the adjoint, while in accord with modern 
usage, differs from the definitions given in some older textbooks. 
A Hermitian matrix is one which is equal to its own adjoint : 



H = H f or Hij = Ha*. 



(67) 



The product of an m X n matrix and an n-component vector is an 
m-component vector, written 
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y = Ax or v/ t - = E AijXj , i = 1,2, ••■ ,m. (68) 

i=i 

If x is an ra-component vector, y an n-component vector, and A an 
m X n matrix, then from (64) and (66), 

(x,Ay) = £ ZxiAjy* = (A f x,y). (69) 

i=i j=i 

A Hermitian form is the scalar product of Hx with x, where H is a 
Hermitian matrix : 

(Hx,x) - E Y.xtHaXi. (70) 

i=i j=i 

Hermitian forms are real-valued, since in view of (67), 

(Hx,x)* = £ £*ff«V - £ txfHjiXi - (Hx,x). (71) 
t=i j=i »"=i j=i 

A Hermitian form is positive definite if 

(Hx,x) > whenever (x,x) ^ 0. (72) 

If the > sign is replaced by ^ , the form is called positive semide finite. 
The product of an m X w matrix A and an n X p matrix B is an 
m X p matrix C whose elements are given by 

Cij = E A ik B ki , i = 1,2, ■■■ ,m; j = 1,2, ■■■ , p. (73) 

A square matrix whose determinant vanishes is called singidar. If the 
determinant does not vanish, the matrix is called nonsingular. The matrix 
of a positive definite Hermitian form is nonsingular. 

The inverse of a nonsingular n X n square matrix A is the n X n 
square matrix A - which satisfies 

A"'A = AA" 1 = 1„ , (74) 

where 1„ is the n X n unit matrix with l's on the main diagonal and 
0's elsewhere. The elements of A~' are given by 

< A -'>* = ai^ (75) 

where GLji is the cof actor of the element A a in the determinant of A. 

If r is an arbitrary m X n matrix, r + r is an n X n Hermitian matrix, 
since 
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(r'r),, = £ r ki *v kj = £ r kj T ki * = (r + r)/, 

A=l A=l (7t>) 

i = 1, 2, ••• ,n; i = 1, 2, ••• , n. 

If A and B are w X n Hermitian matrices and X is a complex parameter, 
then A — XB is called a pencil of matrices. If B is positive definite, the 
pencil is called regular. The characteristic equation of a regular pencil, 
namely 

det (A - XB) = 0, (77) 

always has n real roots \ x , X 2 , • • • , X» , which are called the eigenvalues 
of the pencil. The eigenvalues correspond to eigenvectors z , z , • • ■ , 
z (n) , which satisfy the homogeneous equations 

Az (fc) = X A .Bz a '\ fc= 1, 2, ••• ,n. (78) 

The eigenvectors may be chosen to satisfy 

(B« ( V») =&<>> ( 79 ) 

where 8a is the Kronecker delta. 

The largest eigenvalue X„ of the regular pencil A — XB satisfies 

a (Ax,x) / sn s 

X " = n i?(Bi^)> (80) 

and this maximum is assumed only for eigenvectors of the pencil corre- 
sponding to the eigenvalue X„ . 

APPENDIX B 

Reciprocity Theorem for Time-Harmonic Fields 

We shall prove the reciprocity theorem in a form convenient for use in 
the present paper, following an approach similar to that of Harrington. 5 

Consider a linear, time-invariant medium characterized by the permit- 
tivity tensor e, the permeability tensor y, and the conductivity tensor d. 
All three tensors are assumed to be symmetric, although they may be 
functions of the space coordinates. Let 

<y = d + iccz, Z = «'coy , (81) 

where w is the angular frequency of the time-harmonic fields. 

Consider two sets of electric current densities, J a and J b , which are 
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vector functions of position and which give rise to the fields E°, H a and 
E b , H b respectively. Maxwell's equations are 

V X H" = OjE a + T, V X H" = <yE 6 + J b , 

-V X E" = ZH a , -V X E 6 = ZH 6 . 

From the first and fourth equations, 

-V-(E b X H") = E*-V X H" - H"V X E" 

= E 6 -«yE° + E b -J° + H fl -ZH 6 , 

and from the second and third, 

-V-(E" X H") = E"V X H" - H 6 V X E" 

(84) 

= F-'yE 6 + E u -J 6 + H b ZH". 

Subtracting (83) from (84) and using the symmetry of «y and Z, we 
obtain 

V-(E" X H" - E" X H") = F-J" - E"-r. (85) 

Now integrate over a large spherical volume V bounded by the surface 
S, which contains all sources and matter in its interior. The divergence 
theorem yields 

/ (E b X H° - E a X H") n dS = £ (E a - f - E b -J a ) dV, (86) 

where n is the outward normal to S. The individual fields fall off as 1/r, 
where r is the radius of V, but for large r the leading terms satisfy 

E" = rjH" X n, E" = ^H" X n, (87) 

where i\ is the characteristic impedance of free space. 
Hence for the leading terms, 

rT x [E b X H" - E" XH 6 ]n 

= [(H b X n) X H" - (H a X n) X H b ]-n (88) 

= [n(H a -H b ) - H"(n-H a ) - n(H a -H b ) + H a (n-H b )]-n - 0. 

It follows that if all sources and matter are of finite extent, then 

fE a -fdV= /E b -jW, (89) 

where each integral is taken over the region in which the source currents 
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are different from zero. If the medium is not symmetric, the theorem is 
still true provided that E represents the field produced by J in the 
"transposed" medium; but this generalization is not very useful in the 
present context. 

Now let J° correspond to an electric dipole of unit moment in the di- 
rection u" at the point P a , and let J correspond to an electric dipole of 
unit moment in the direction u at P b • Equation (89) takes the form 

u 6 -E"(P 6 ) = u a -E''(P„), (90) 

where the left side represents the components of electric field due to 
source A at the location and in the direction of source B, and the right 
side represents the component due to source B at the location and in the 
direction of source A . This is the desired reciprocity theorem. 



appendix c 

Sandberg's Lemma 

We reproduce Sandberg's proof 3 of the following result. 

Lemma. If A and B respectively are n X m and m X n matrices, then 

det (1„ + AB) = det (l m + BA). 

Proof. First consider the case in which A and B are square p X p 
matrices. Then, if A is nonsingular, 

det [l p + AB] = det [A -1 (l p + AB)A] 

(91) 
= det [l p + BA]. 

If A is singular, it has a zero characteristic root, and hence there exists a 
positive number X u such that A + M p is nonsingular for all real X satis- 
fying < | X | < Ao . Thus when < | X | < X n , 

det [1„ + (A + X1 P )B] = det [1, + B(A + XI,)]. (92) 

Both sides of (92) are polynomials in X of degree at most p. Furthermore 
these polynomials must be identical since they agree throughout the 
real interval (0,Xo). Therefore (92) is valid when X = 0. 

Consider now the case in which A and B are not square. Let 
p = m + n, 



m n n m 

A |_0 0_K |_0 Ojn' 



(93) 
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and let the symbol -j- denote a direct sum of matrices. Observe that 

det [l p + AB] = det [(1„ + AB) -j- l m ] = det [1„ + AB], 

(94) 
det \l p + BA] = det [(1„, + BA) + 1„] = det [l m + BA], 

which proves the lemma. 
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